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LONG-TERM  GOALS 

Develop  a  simple  parameterization  for  the  direetional  wave  speetrum  in  the  vicinity  of  a  hurricane. 

OBJECTIVES 

Develop  and/or  modify  the  real-time  operating  system  and  analysis  techniques  and  programs  of  the 
NASA  Seanning  Radar  Altimeter  (SRA)  to  proeess  the  SRA  wave  topography  data  into  direetional 
wave  spectra  during  hurricane  flights.  Upload  the  speetra  and  the  topography  onto  a  web  site 
immediately  post-flight  to  make  them  available  to  ONR  investigators. 

APPROACH 

The  SRA  has  a  long  heritage  in  measuring  the  energetie  portion  of  the  sea  surface  direetional  wave 
speetrum  (Walsh  et  al.  1985;  1989;  1996,  2002;  Wright  et  al.  2001).  The  wave  speetra  have  been  used 
to  assess  the  output  of  the  WaveWateh  III  numerieal  wave  model  (Moon  et  al.,  2003;  Fan  et  al.,  2007). 
To  obtain  the  direetional  wave  speetrum,  the  energy  in  the  eneounter  speetrum  generated  from  the 
SRA  wave  topography  must  be  doubled  everywhere,  the  artifact  lobes  deleted,  and  the  real  lobes 
Doppler-correeted.  Identifying  the  artifaet  lobes  for  deletion  and  partitioning  the  real  speetral  lobes 
into  the  various  wave  eomponents  has  been  a  slow  and  labor-intensive  proeess  and  efforts  have  been 
made  toward  eliminating  operator  deeisions  in  that  proeess  to  enable  unattended  operation. 

WORK  COMPLETED 

After  flying  into  hurrieanes  aboard  a  NOAA  aireraft  for  eight  seasons  (1998-2005)  the  NASA  SRA 
experienced  significant  hardware  problems  and  was  deeommissioned.  NOAA  awarded  a  eontract  to 
ProSensing  (prosensing. eom)  to  build  an  operational  SRA.  It  should  be  shipped  to  the  NOAA  Aireraft 
Operations  Center  in  Tampa  by  15  October  2007.  The  NASA-developed  software  and  analysis 
teehniques  have  been  tailored  to  the  new  NOAA  SRA  hardware  and  system  eharacteristics  to  enable 
data  proeessing  and  transmitting  final  data  products  to  the  National  Hurricane  Center  (NHC)  during 
the  flights  to  provide  an  enduring  legaey  to  ONR  investigators  and  the  nation. 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 
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RESULTS 


Figure  1  shows  the  output  of  the  initial  SRA  display  software  which  has  been  installed  and  tested  on  an 
NHC  computer  using  data  acquired  during  a  flight  into  Hurricane  Humberto  on  24  September  2001. 
Each  time  the  program  is  run,  it  reads  in  all  the  SRA  data  files  sent  from  the  aircraft  and  accumulated 
in  a  directory  in  the  NHC  computer,  finds  the  maximum  SWH  and  automatically  generates  a  plan  view 
of  the  aircraft  track  through  the  hurricane,  color  coding  the  wave  height  red  at  the  16',  20',  24',  28',  32', 
or  36'  level,  depending  on  the  maximum  level  the  wave  height  exceeds.  So  unless  the  waves  are  lower 
than  16',  the  display  will  always  contain  all  the  colors  (green,  black,  blue,  red).  Future  programs  will 
display  entire  directional  spectra. 
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Figure  1.  Display  for  SRA  software  installed  and  tested  on  National  Hurricane  Center  computer. 


The  2006  annual  report  (Document  Number  N0001406IP20037)  included  an  analysis  focused  on 
quantifying  distortions  in  the  measured  wave  spectra  produced  by  variations  in  the  power 
backscattered  by  the  sea  surface.  The  antenna  of  the  new  NOAA  SRA  forms  its  beams  digitally 
instead  of  scaiming  a  lens-generated  beam  mechanically  as  the  NASA  system  did,  but  both  two-way 
beamwidths  are  1°  and  the  effective  swath  is  between  22°  off-nadir  to  the  left  of  the  aircraft  track, 
through  nadir,  to  22°  off-nadir  to  the  right.  The  solid,  parallel  lines  in  Figure  2  (from  the  2006  report) 
show  the  0.9  m  half-power  width  bounds  of  the  SRA  pulse  approaching  the  sea  surface  from  1500  m 
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altitude.  The  horizontal  sinusoidal  eurve  represents  an  oeean  wave  of  50  m  length  and  1.67  m  erest-to- 
trough  height.  The  solid  line  orthogonal  to  the  transmitted  pulse  indieates  the  antenna  boresight  and 
the  dashed  and  dotted  lines  parallel  to  it  indieate  the  half-power  and  0.1  power  limits  of  the  1°  two- 
way  antenna  pattern,  modeled  as  Gaussian. 


Figure  2.  Nadir  and  off-nadir  geometry  for  the  SRA  pulse  transmitted 
from  1500  m  height  interacting  with  an  ocean  wave  of  50  m  length. 


It  is  apparent  that  the  eentroid  of  the  baekseattered  power  for  the  pulse  approaehing  the  nadir  point  in 
the  top  panel  will  determine  a  range  greater  than  that  to  the  wave  erest  at  the  boresight  beeause  the 
half-power  width  of  the  antenna  pattern  extends  to  mean  sea  level.  The  reverse  would  be  true  for  a 
trough  at  the  boresight.  This  spatial  filtering  by  the  antenna  footprint  reduees  the  measured  wave 
height  at  nadir. 

In  the  bottom  panel  of  Figure  2  the  antenna  boresight  is  20°  off-nadir.  The  wave  erest  at  the  boresight 
makes  the  loeal  incidenee  angle  20°  there,  but  places  the  maximum  ±6°  slopes  of  the  50  m  length  wave 
at  the  antenna  half-power  points.  The  local  incidence  angle  is  26.5°  at  the  off-nadir  antenna  half¬ 
power  limit  and  13.5°  at  the  near-nadir  limit.  The  smaller  local  incidence  angle  on  the  near-nadir  side 
of  the  wave  crest  will  increase  the  baekseattered  power  and  shift  its  centroid  to  shorter  range. 

Ascribing  the  shorter  range  to  the  antenna  boresight  angle  will  increase  the  apparent  height  of  the  wave 
crest.  If  a  wave  trough  were  at  the  antenna  boresight,  the  local  incidence  angle  would  be  15.5°  at  the 
off-nadir  half-power  antenna  limit  and  25.5°  at  the  near-nadir  limit.  The  lower  local  incidence  angle 
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on  the  off-nadir  side  of  the  trough  would  inerease  the  baekseattered  power  and  shift  the  power  eentroid 
to  longer  range,  erroneously  indieating  a  deeper  wave  trough  at  the  antenna  boresight. 

To  be  able  to  model  and  eorreet  these  distortions  it  is  neeessary  to  know  the  variation  of  the 
baekseattered  power  with  ineidence  angle,  whieh  is  a  function  of  the  sea  surface  mean  square  slope 
(mss),  principally  determined  by  the  small  scale  roughness.  The  Plant  (1982)  theoretical  maximum 
limit  on  mss  at  high  wind  speed  is  0.080.  The  baekseattered  power  is  generally  assumed  to  fall  off  in  a 
Gaussian  fashion  near  nadir  with  the  falloff  rate  being  determined  by  the  reciprocal  of  the  mss  (Walsh 
et  al.  2008). 

When  the  log  of  the  baekseattered  power  is  plotted  versus  the  square  of  the  incidence  angle  tangent, 
the  Gaussian  assumption  produces  a  straight  line.  But  Figure  3  shows  SRA  data  indicating  that  its 
measurements  extend  far  enough  off  nadir  that  the  strictly  Gaussian  assumption  is  not  valid.  The  data 
are  better  fit  with  curves  of  the  form 


lnP(S)  =  -AS^  +BS^  (1) 

where  S  is  the  incidence  angle  tangent.  In  the  Gaussian  assumption,  B  =  0.  The  B  term  moderates  the 
falloff  of  the  baekseattered  power  and  reduces  the  distortion  in  the  spectral  measurements. 
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SQUARE  OF  LOCAL  INCIDENCE  ANGLE  TANGENT 

Figure  3.  Backscattered  power  variation  for  the  SRA  scan  plane  within  ±7.5°  of  the  upwind- 
downwind  direction  with  positive  angles  indicating  the  downwind  direction.  Circles  indicate  data 
for  which  the  right  side  of  the  SRA  scan  plane  was  oriented  downwind  and  +  symbols  indicate  the 
data  for  which  the  right  side  of  the  scan  plane  was  oriented  upwind.  The  data  for  8  m  s'^  were 
shifted  down  by  10  dB  and  the  data  for  2  m  s'^  were  shifted  down  by  15  dB  to  prevent  overlap.  The 
second-order  curves  on  the  right  side  were  least-squares  fitted  to  the  circles  and  the  curves  on  the 

left  side  were  least-squares  fitted  to  the  +  signs. 


There  were  several  problems  in  trying  to  use  SRA  hurrieane  data  to  investigate  the  baekseattered 
power  variation  with  ineidenee  angle.  At  the  generally  high  altitudes  of  the  hurrieane  flights  the  SRA 
signal  level  was  low,  the  footprint  was  large,  and  data  were  not  available  over  a  wide  variety  of 
azimuths  to  do  a  eomprehensive  analysis.  The  SRA  data  used  were  aequired  at  low  altitude  in  the 
Southern  Oeean  with  eomplete  azimuthal  eoverage  over  a  wide  range  of  sea  states  (up  to  9  m  wave 
height)  and  wind  speeds,  though  less  than  hurrieane  foree  (Walsh  et  al.  2008).  Figure  4  shows  the 
measurement  geometry.  Data  were  aequired  with  the  aireraft  in  a  7°  roll  attitude  while  it  turned 
azimuthally  through  810°.  The  small  antenna  footprint  (4  m)  allowed  eomputation  of  the  along-  and 
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cross-track  components  of  sea  surfaee  slope  at  eaeh  measurement  point  to  eompute  loeal  ineidenee 
angles  for  eaeh  baekseattered  power  value. 


Figure  4.  Measurement  geometry  for  data  acquired  at 
250  m  altitude  while  the  aircraft  was  in  a -  7"  roll  attitude. 
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Figure  5  shows  a  representative  data  segment.  Aequirmg  data  at  all  azimuths  relative  to  the  wind 
direetion  required  spending  relatively  little  time  in  any  partieular  azimuth  interval.  The  slopes 
assoeiated  with  the  few  dominant  waves  observed  in  a  particular  azimuth  interval  would  not  be 
statistically  representative  of  the  overall  slope  distribution  of  the  sea  state  and  would  not  even  be  the 
same  at  all  off-nadir  incidence  angles  in  any  particular  azimuth  interval.  The  result  would  be  a 
significant  misbinning  of  the  rapidly  fluctuating  backscattered  power  data  with  respect  to  the  actual 
incidence  angles  and  a  corruption  of  the  resulting  distributions.  The  high  quality  of  the  surface  slope 
data  displayed  in  Figure  5  suggest  that  sea  surface  slopes  on  scales  too  small  to  be  resolved  by  the 
SRA  wave  topography  measurements  would  be  small  enough  that  their  distribution  could  be 
considered  statistically  identical  at  all  incidence  angle  and  azimuth  intervals. 
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Figure  5.  SRA  wave  topography  (left  panel)  and  backscattered  power  with  the  average  cross-track 
power  variation  removed  (second  panel).  The  white  line  indicates  the  nadir  point.  The  third  and 
fourth  panels  are  the  along-track  and  cross-track  sea  surface  slopes,  respectively. 


The  circles  in  Figure  6  show  the  quadratic  (B)  coefficient  plotted  against  the  linear  (A)  coefficient 
(equation  1)  for  all  the  second-order  fits  to  the  SRA  data  such  as  shown  by  the  curves  in  Figure  3.  The 
straight  line  in  Figure  6  was  least-squares  fitted  to  the  logarithms  of  the  B  and  A  coefficients  for  all  the 
data.  Its  equation  is  the  power  law 


5  =  0.4182^^-^^^  ^2) 

The  straight  line  with  dots  was  least-squares  fitted  to  the  logarithms  of  the  B  and  A  coefficients  for  all 
the  data  except  the  lightest  wind-speed  day  (60<y4<100),  resulting  in 


B  =  0.5616 

The  dashed  curve  is  actually  a  straight  line  (appearing  curved  due  to  the  log-log  axes  in  Figure  6)  that 
was  least-squares  fitted  to  the  B  and  A  coefficients  for  all  the  data  except  the  lightest  wind-speed  day. 
Its  equation  is 


5  =  -9.84 +  2.066^,  (4) 

The  dashed  and  solid  curves  indicate  a  very  similar  dependence  of  5  on  ^4  for  values  of  A  less  than  40. 
If  the  cluster  of  data  for  the  2  m  s"'  wind  speed  day  (60<y4<100)  were  included  in  the  linear 
dependence  curve  fit  (4),  the  fit  to  the  main  body  of  data  would  have  degraded  significantly.  The 
equation  of  the  dotted  line,  least-squares  fitted  to  the  log  of  the  three  Ku-band  values  (+  symbols)  from 
the  Voronovich  and  Zavorotny  (2001)  theoretical  analysis  is 


5  =  0.4194^^'^^^ 

This  analysis  indicates  that  the  non-Gaussian  characteristic  of  the  scattering  increases  as  the  mss 
increases.  Figure  6  suggests  the  possibility  that  there  may  be  a  global  model  to  characterize  Ka-band 
radar  backscatter  from  the  sea  surface  within  25°  of  nadir  and  that  the  quadratic  {B)  coefficient  is 
uniquely  related  to  the  linear  (A)  coefficient.  Except  for  the  lightest  wind  day,  there  was  considerable 
overlap  in  the  A  and  B  coefficients  from  day  to  day  because  they  varied  significantly  azimuthally.  But 
the  scatter  for  any  given  day  spread  the  coefficients  along  the  general  power  law  representing  the 
entire  dataset.  This  is  important  because  it  suggests  that  for  a  particular  value  of  the  A  coefficient,  the 
same  B  coefficient  will  be  associated  with  it  whether  the  A  coefficient  corresponded  to  the  crosswind 
direction  at  a  high  wind  speed,  or  to  the  upwind  direction  at  a  lower  wind  speed,  or  to  an  intermediate 
wind  speed  at  an  intermediate  direction  relative  to  the  wind. 

The  Voronovich  and  Zavorotny  (2001)  Ku-band  theoretical  analysis,  which  used  only  Gaussian  slope 
statistics  but  included  diffraction,  produced  almost  the  same  power  law  relationship  between  A  and  B 
as  the  SRA  Ka-band  observations.  Further  investigation  is  needed  to  determine  whether  the  offset 
between  the  Ku-band  curve  fit  to  their  theoretical  results  and  the  Ka-band  data  is  due  to  a  contribution 
of  non-Gaussian  slope  statistics,  or  to  the  frequency  difference,  or  to  problems  in  the  theory  and/or 
measurements.  But  what  it  suggests  is  that  whether  the  non-Gaussian  nature  of  the  scattering  is  caused 
by  non-Gaussian  slope  statistics,  or  Gaussian  slope  statistics  and  diffraction,  or  some  combination  of 
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the  two,  the  linear  coefficient,  A,  is  sufficient  to  determine  the  characteristics  of  the  backscatter  since 
the  quadratic  coefficient,  B,  can  be  related  to  it. 

It  should  be  emphasized  that  the  relationship  between  A  and  B  was  not  obtained  using  off-nadir  angles, 
but  local  incidence  angles  which  took  into  account  the  slopes  of  waves  whose  wavelengths  were 
longer  than  about  1 6  to  25  m.  Equation  (2)  resulting  from  this  analysis  can  be  used  in  the  correction  of 
SRA  hurricane  spectra  but  it  may  have  more  general  application.  Computing  local  incidence  angles  is 
not  possible  for  satellite  remote  sensors.  But  their  measurements,  binned  by  off-nadir  angle,  are 
generally  averaged  over  areas  large  enough  that  long-wave  slopes  can  be  assumed  to  be  statistically 
represented  at  all  incidence  angles.  This  simple  model  might  prove  very  useful  for  predicting  the 
performance  of  remote  sensors  that  operate  within  25°  of  nadir. 


Figure  6.  Quadratic  (B)  coefficient  versus  the  linear  (A)  coefficient  (1)  for  all  the 
second  order  fits  to  the  SRA  data  of  Figure  3.  The  straight  line  is  a  power  law  least 
-squares  fitted  to  all  the  data.  The  straight  line  with  dots  is  power  law  least-squares  fitted 
to  the  data  for  which  A  <  40.  The  dashed  curve  is  a  straight  line  least-squares  fitted  to 
the  data  for  which  A  <  40.  The  dotted  straight  line  is  a  power  law  least-squares  fitted  to 
the  three  values  (+)  from  the  Voronovich  and  Zavorotny  (2001)  Ku-band  theoretical  analysis. 
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IMPACT/APPLICATIONS 


The  SRA  provided  the  first  comprehensive,  quantified  measurements  of  the  directional  wave  spectrum 
spatial  variation  in  the  vicinity  of  hurricanes.  The  data  impact  all  assessments  of  air/sea  interaction  in 
a  hurricane  environment  and  serve  to  validate  wave  models.  The  ability  to  examine  the  structure  of 
individual  waves  and  wave  groups  is  important  for  assessing  viability  of  marine  structures. 

RELATED  PROJECTS 

All  hurricane  components  of  ONR  CBLAST. 
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